Numerical study of Cosmic Ray Diffusion in MHD turbulence

We study diffusion of Cosmic Rays (CRs) in turbulent magnetic fields using test particle simulations. Electromagnetic fields are produced in direct numerical MHD simulations of turbulence and used as an input for particle tracing, particle feedback on turbulence being ignored. Statistical transport coefficients from the test particle runs are compared with earlier analytical predictions. We find qualitative correspondence between them in various aspects of CR diffusion. In the incompressible case, that we consider in this paper, the dominant scattering mechanism occurs to be the non-resonant mirror interactions with the slow-mode perturbations. Perpendicular transport roughly agrees with being produced by magnetic field wandering.Comment: 8 pages, 8 figures, ApJ sub

Comparison of spectral slopes of magnetohydrodynamic and hydrodynamic turbulence and measurements of alignment effects

We performed a series of high-resolution (up to 1024^3) direct numerical simulations of hydro and MHD strong turbulence. We found that for simulations with normal viscosity the slopes for spectra of MHD are similar, although slightly more shallower than for hydro simulations. However, for simulations with hyperviscosity the slopes were very different, for instance, the slopes for hydro simulations showed pronounced and well-defined bottleneck effect, while the MHD slopes were relatively much less affected. We believe that this is indicative of MHD strong turbulence being less local than Kolmogorov turbulence. This calls for revision of MHD strong turbulence models that assume local "as-in-hydro case" cascading. Nonlocality of MHD turbulence casts doubt on numerical determination of the slopes with currently available (512^3--1024^3) numerical resolutions, including simulations with normal viscosity. We also measure various so-called alignment effects and discuss their influence on the turbulent cascade.Comment: 10 pages, 6 figures, extended version, ApJ accepte

Density scaling and anisotropy in supersonic MHD turbulence

We study the statistics of density for supersonic turbulence in a medium with magnetic pressure larger than the gaseous pressure. This study is motivated by molecular cloud research. Our simulations exhibit clumpy density structures, which contrast increases with the Mach number. At 10 Machs densities of some clumps are three orders of magnitude higher than the mean density. These clumps give rise to flat and approximately isotropic density spectrum corresponding to the random distribution of clumps in space. We claim that the clumps originate from our random, isotropic turbulence driving. When the contribution from those clumps is suppressed by studying logarithm of density, the density statistics exhibit scale-dependent anisotropy consistent with the models where density structures arise from shearing by Alfv\'en waves. It is noteworthy that originally such models were advocated for the case of low-Mach, nearly incompressible turbulence.Comment: 4 pages, 4 figures, using emulateap

Growth of Magnetic Fields Induced by Turbulent Motions

We present numerical simulations of driven magnetohydrodynamic (MHD) turbulence with weak/moderate imposed magnetic fields. The main goal is to clarify dynamics of magnetic field growth. We also investigate the effects of the imposed magnetic fields on the MHD turbulence, including, as a limit, the case of zero external field. Our findings are as follows. First, when we start off simulations with weak mean magnetic field only (or with small scale random field with zero imposed field), we observe that there is a stage at which magnetic energy density grows linearly with time. Runs with different numerical resolutions and/or different simulation parameters show consistent results for the growth rate at the linear stage. Second, we find that, when the strength of the external field increases, the equilibrium kinetic energy density drops by roughly the product of the rms velocity and the strength of the external field. The equilibrium magnetic energy density rises by roughly the same amount. Third, when the external magnetic field is not very strong (say, less than ~0.2 times the rms velocity when measured in the units of Alfven speed), the turbulence at large scales remains statistically isotropic, i.e. there is no apparent global anisotropy of order B_0/v. We discuss implications of our results on astrophysical fluids.Comment: 16 pages, 18 figures; ApJ, accepte

Cosmic Ray transport through gyroresonance instability in compressible turbulence

We study the nonlinear growth of kinetic gyroresonance instability of cosmic rays (CRs) induced by large scale compressible turbulence. This feedback of cosmic rays on turbulence was shown to induce an important scattering mechanism in addition to direct interaction with the compressible turbulence. The linear growth is bound to saturate due to the wave-particle interactions. By balancing increase of CR anisotropy via the large scale compression and its decrease via the wave-particle scattering, we find the steady state solutions. The nonlinear suppression due to the wave-particle scattering limit the energy range of CRs that can excite the instabilities and be scattered by the induced slab waves. The direct interaction with large scale compressible modes still appears to be the dominant mechanism for isotropization of high energy cosmic rays (> 100 GeV).Comment: 11 pages, 6 figures, 1 table, typos corrected, ApJ in pres

The Spectral Slope and Kolmogorov Constant of MHD turbulence

The spectral slope of strong MHD turbulence has recently been a matter of controversy. While Goldreich-Sridhar model (1995) predicts Kolmogorov's -5/3 slope of turbulence, shallower slopes were often reported by numerical studies. We argue that earlier numerics was affected by driving due to a diffuse locality of energy transfer in MHD case. Our highest-resolution simulation (3072^2x1024) has been able to reach the asymptotic -5/3 regime of the energy slope. Additionally, we found that so-called dynamic alignment, proposed in the model with -3/2 slope, saturates and therefore can not affect asymptotic slope. The observation of the asymptotic regime allowed us to measure Kolmogorov constant C_KA=3.2+-0.2 for purely Alfv\'enic turbulence and C_K=4.1+-0.3 for full MHD turbulence. These values are much higher than the hydrodynamic value of 1.64. The larger value of Kolmogorov constant is an indication of a fairly inefficient energy transfer and, as we show in this Letter, is in theoretical agreement with our observation of diffuse locality. We also explain what has been missing in numerical studies that reported shallower slopes.Comment: 5 pages 3 figure